ggplot2 Kunduz1401

library(rgdal)
library(ggplot2)
library(geojsonio)
library(dplyr)
library(broom)
library("arules")
library("arulesViz")
library(mapproj)
library("ggrepel")

afghanistan.json <- geojson_read("jsons/district398.json", what="sp")
k_event_df <- read.csv("data/k_event.csv", header = TRUE)
  
k_event_df

plot(afghanistan.json)

afghanistan.json@data

afghanistan_fortified <- fortify(afghanistan.json, region = "DISTID")

Too few points in geometry component at or near point 60.928091700000003 33.497045180000001SpP is invalid
Invalid objects found; consider using set_RGEOS_CheckValidity(2L)

map <- ggplot(data = afghanistan_fortified, aes(x = long, y = lat, group = group))
map + geom_path()

map + 
  geom_polygon(aes(fill = id)) +
  coord_fixed(1.3) +
  guides(fill = FALSE)

afghanistan_fortified %>% distinct(id) %>% write.csv("data/districts.csv", row.names= FALSE)

cities <- data.frame("name" = c("Kabul", "Kandahar", "Herat", "Mazar-i-Sharif"), "long" = c(69.160652, 65.7372, 62.2040, 67.1180), "lat" = c(34.543896, 31.6289, 34.3529, 36.6926))

add_k_to_map <- function(k_days) {
  k_days_df <- k_event_df[which(k_event_df$visit_day - k_event_df$impact_day == k_days),]
  print(length(k_days_df))
  map_districts <- read.csv("data/districts.csv")
  map_districts$visits = 0
  for (i in 1:nrow(k_days_df)) {
    map_districts$visits[map_districts$id == k_days_df[i, 'destination_district']] <- k_days_df[i, 'visits']
  }
  
  # plot the shaded map_districts
  afghanistan_fortified_days <- merge(afghanistan_fortified, map_districts, by="id")
  map <- ggplot(data = afghanistan_fortified_days, aes(x = long, y = lat, group = group))
  
  # create the legend, title, and caption for the map
  subt <- paste("Number of impacted people who were in another district", k_days, "day(s) after the day of violence.")
  map <- map + geom_polygon(aes(fill=visits), color='grey', size=.1) + scale_fill_gradient(high = "#e34a33", low = "#ffffff", limits=c(0, 3698), guide = "colorbar")  + coord_fixed(1.3)
  map <- map + labs(title = "Displacement following the Battle of Kunduz", subtitle=subt) + theme(plot.title = element_text(size=20, face="bold"))
  
  
  # make the Kunduz district (id: 1401) a separate color from the rest
  map <- map + geom_polygon(data = afghanistan_fortified_days[which(afghanistan_fortified_days$id == 1401),], aes(fill=visits), color='black', size=.1, fill = "blue")
  
  # add the cities: Kabul, Kandahar, Herat, Mazar-i-Sharif
  map <- map + geom_point(data=cities, aes(x = long, y = lat, group = name),col="black", size=1) + geom_text_repel(data=cities, segment.size =.2, aes(x = long, y = lat, label= name, group=name), fontface= "bold", size=4, color = "black", nudge_x = c(3, 2, -2, 1), nudge_y = c(0.25, -1, 0.5, .8))
  
  map
}
  
add_k_to_map(1)

[1] 6

add_k_to_map(5)

[1] 6

add_k_to_map(10)

[1] 6

add_k_to_map(15)

[1] 6

add_k_to_map(30)

[1] 6

add_k_to_map(45)

[1] 6

add_k_to_map(60)

[1] 6

add_k_to_map(75)

[1] 6

add_k_to_map(90)

[1] 6

add_k_to_map(105)

[1] 6

add_k_to_map(120)

[1] 6

IGNORE EVERYTHING BELOW THIS LINE OF CODE!!!!

map <- ggplot(data = afghanistan_fortified, aes(x = long, y = lat, group = group))

map + geom_polygon(aes(fill=visits), color='black', size=.1) + scale_fill_gradient(high = "#e34a33", low = "#ffffff", limits=c(0, 3386), guide = "colorbar")  + coord_fixed(1.3)

k_days_df <- k_event_df[which(k_event_df$visit_day - k_event_df$impact_day == 15),]
for (i in 1:nrow(k_days_df)) {
  #print(k15_df[i, 'visits'])
  afghanistan.json@data$visits[afghanistan.json@data$DISTID == k_days_df[i, 'destination_district']] <- k_days_df[i, 'visits']
}
afghanistan.json@data
#plot(afghanistan, col=colors, bg="#A6CAE0", cex.main=3)

afghanistan.json_fortified <- tidy(afghanistan.json, region = "visits")
afghanistan.json_fortified

ggplot() +
  geom_polygon(data = afghanistan.json_fortified, aes( x = long, y = lat, group = group), fill="white", color="black") +
  theme_void() +
  coord_map()